NYO-1480-93 


Courant  Institute  of 
Mathematical  Sciences 

AEG  Computing  and  Applied  Mathematics  Center 


Non-Empirical  SCF-MO  Calculations 
on  Transition  Metal  Complexes 

Harold  Basch,  C.  Hollister, 
and  Jules  W.  Moskowitz 


AEC  Research  and  Development  Report 

Mathematics  and  Computing 
July  1968 


New  York  University 


NtlVv'  YOi<J<  UNIVGi<.:>iTY 
COURANT  INSTITUTE  -  LltiRARY 
251  Mercer  St.    New  York,  NY.  10012 


UNCLASSIFIED 


AEC  Computing  and  Applied  Mathematics  Center 
Courant  Institute  of  Mathematical  Sciences 
New  York  University 


Mathematics  and  Computing  NYO-l480-93 

NON-EMPIRICAL  SGF-MO  CALCULATIONS  ON  TRANSITION  METAL  COMPLEXES 

Harold  Basch,  C.  Holllster  and  Jules  W.  Moskowltz 


Contract  No.  AT(30-l)-l48o 


UNCLASSIFIED 


In  this  paper  we  shall  demonstrate  the  feasibility  of 
employing  Gaussian  fui ctlons  as  basis  orbltals  in  molecular 
orbital  computations  on  transition  metal  complexes.   This 
calculation  represents  the  first  treatment  of  a  transitlor, 
metal  complex  ion  involving  all  electrons  and  no  approximations 
within  the  Roothaan  SCF-MO  formalism.   The  qualitative  features 
of  the  calculation  can  be  expected  to  have  wide  significance 
both  for  the  understanding  of  the  electronic  structure  of 
transition  metal  complexes  and  for  the  application  of  MO 
theory  to  large  systems  of  this  type. 

We  will  report  the  results  of  calculations  on  the  (hypo- 

o  _ 

thetical)  square  planar  NlFu    ion  (the  reasons  for  this 
particular  choice  will  be  explained  later)  and  discuss  those 
irteresting  ground  and  excited  state  molecular  electronic 
properties  which  can  be  extracted  from  a  computed  wave  functlor. 

It  is  well  known  that  individual  Gaussians  are  much 
poorer  representations  of  atomic  orbltals  (aO's)  than  are 
single  Slater  type  orbltals  (STO's),  especially  in  two  very 
Important  areas:  near  the  nucleus,  and  in  the  tail  region  of 
the  AO.   These  deficiencies  are  very  serious  if  ore  warits 
to  use  Gaussians  as  basis  functions  for  molecular  calculations. 
Of  course,  a  single  STO  is  also  kr  owr.  to  reproduce  rather 
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poorly  the  tail  regioii  of  an  AO,  and  It  takes  a  linear  combina- 
tion of  STO's  to  span  adequately  the  whole  space  of  an  AO.   This 
suggests  the  method  whereby  one  can  correct  the  deficiencies 
of  the  individual  Gaussian  type  orbitals  (GTO's);  use  many 
of  them  in  linear  combinations  to  form  Gaussian  type  functions 
(GTF's)  and  then  use  the  GTP's  as  the  basis  functions  in  the 
MO  calculation.   The  proper  linear  combinations  needed  to  form 
the  GTF's  are  easily  obtained  either  by  directly  fitting  the 
Hartree-Fock  AG's  with  GTO's,  or  doing  an  atomic  analytic  SCF 
calculation  in  a  GTO  basis.   We  find  that  on  the  average  only 
2^-3  GTO's  per  STO  are  required  for  quantitatively  comparable 
results . 

This  is  Illustrated  in  Figure  I  which  compares  the  3d  AO 

p-J-        Q 

of  Ni    (3d   configuration)  in  the  single  Slater,  Hartree-Fock, 
and  Gaussian  approximations.   Note  that  the  position  of  the 
maximum  for  the  Hartree-Fock  and  GTF  coincide  and  that  they 
tail  off  pretty  much  together,  in  contrast  to  the  single  Slater 
function.   The  Gaussian  curve  is  a  linear  combination  of  5  GTO's 

and  is  identical,  on  the  scale  of  this  figure,  with  the  2  term 

f2) 
Slater  representation  computed  by  J.W.  Richardson,^  '' 

This  ratio  of  2-3  Gaussians  per  Slater  means  that  we  must 

use  a  basis  set  2-3  times  larger  for  comparable  results  and 
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evaluate  many  more  integrals.   Fortunately,  the  speed  with 
which  the  necessary  two-electron  integrals  over  GTO's  can  be 
computed  more  than  compensates  for  the  handicap  of  the  larger 
basis.   This  very  favorable  speed  is  due  entirely  to  the  fact 
that  all  the  multicenter  integrals  over  GTO's  can  be  evaluated 
using  simple  analytic  formulas  which  are  easily  coded  In 
FORTRAN  language,  without  recourse  to  complex  numerical  integra- 
tion techniques,  numerous  slowly  convergent  expansions,  or 

(3) 

sophisticated  programming  structure.^  ^   This  is  the  fundamental 

reason  for  using  Gausslans. 

To  summarize,  the  MO  wave  function  is  an  anti-symmetrized 
product  of  MO's  which  are  expressed  as  linear  combinations  of 
GTF's.   The  coefficients  of  the  GTF's  in  the  MO's  are  deter- 
mined variationally  via  the  well-known  Roothaan  Hartree-Fock 
equations.^  ■'   The  GTF's  are  themselves  linear  combinations  of 
GTO's  but  with  the  coefficients  in  these  linear  combinations 
fixed  either  directly  or  indirectly  by  previous  atomic  SCF 
calculation.   Finally,  and  this  is  to  be  emphasized,  no  loss 
in  accuracy  in  the  computed  MO's  is  necessary  using  Gausslans 
as  basis  functions  provided  that  they  are  used  properly,  in 
the  sense  that  a  sufficient  number  of  them  are  used  to  span 
adequately  the  function  space. 


-3- 


The  transition  metal  Ion  complex  whose  ground  state  wave 
function  we  have  computed  is  the  66  electron  NlFh    ion  in  a 

square  planar  geometry.   Actually,  this  ion  has  not  been  ob- 

2- 
served  in  a  planar  configuration.   However,  the  planar  PdCl^ 

(^) 
complex  ion  is  well  known  and  has  been  very  extensively  studied.^  ' 

Making  reasonable  extrapolations  from  the  known  shifts  in  behavior 

in  going  from  Ni(ll)  to  Pd(ll)  and  F  to  CI   it  was  hoped  to 

obtain  reasonable  good  estimates  of  properties  to  compare  with 

2" 
those  computed  here  for  NiFi^   .   Unfortunately,  this  expecta- 
tion was  not  realized  in  a  quantitative  manner. 

Outside  of  a  specific  Interest  in  the  square  planar  situation, 
the  reasons  this  particular  ion  was  chosen  for  our  Initial  study 
are  three-fold.   Square  planar  Ni(ll)  complexes  are  predicted 
to  be  closed  shell  in  the  ground  state.   This  avoids  the  con- 
ceptual and  practical  difficulties  that  are  known  to  arise  in 
rigorously  applying  the  SCF  method  to  open  shell  electronic 
configurations.   Further,  since  the  Dj,,  symmetry  has  only  2 
doubly-degenerate  representations,  it  is  likely  that  the  ex- 
cited electronic  states  calculated  in  the  virtual  orbital 
approximation  would  not  involve  the  separation  of  multiplet 
states  arising  from  electrons  distributed  in  two  partially 
occupied,  multiply  degenerate  MO's.   Finally,  square  planar 
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2~ 
NlP^i    is  planned  to  be  part  of  a  series  of  calculations  on 

p  — 
Nl(ll)-F  complexes  which  will  include  NiF^    (tetrahedral ) , 

NlFp  (linear),  and  NlFg    (octahedral). 

The  coordinate  system  used  In  the  calculation  has  the 

nickel  atom  at  the  origin  and  the  four  fluorine  atoms  in  the 

o 
xy -plane .   A  Nl-F  bond  length  of  2.00  A  was  assumed. 

2" 
The  basis  set  for  the  square  planar  NlPj,    calculation  is 

given  in  Table  I  which  shows  the  distribution  of  both  the  GTF's 
and  GTO's  in  the  basis  set  with  respect  to  the  AO's  on  Ni  and  F. 
Notice  that  some  of  the  AO's  are  represented  by  two  basis  func- 
tions.  Especially  Important  are  the  3d  on  Nl  and  2p  of  F. 
This  double  zeta  representation  allows  sufficient  freedom  so 
that  the  free  ion  AO's  may  adjust  to  the  molecular  environment 

in  the  course  of  the  SCF  calculation.   For  example,  the  five 

2+ 
term  Gaussian  representation  of  the  Ni   3d  AO  was  split  up 

into  a  long  range  and  a  short  range  component.   This  allows 
the  3d  AO  either  to  expand  or  contract;  i.e.  it  allows  the 
energetics  of  the  calculation  to  determine  the  spatial  exten- 
sion of  the  3d  AO  in  the  complex.   An  analogous  argument  applies 
to  the  F~  2p  which  is  represented  here  by  2p  and  2p'  basis 
functions. 

The  resultant  ground  state  ordering  and  orbital  energies 
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of  the  computed  MO's  are  shown  in  Figure  II.   We  would  hasten 
to  point  out  that  all  previous  notions  of  the  MO  orderings  in 
the  diamagnetic  square  planar  complexes  of  necessity  arose 
from  absorption  spectra  measurements,  whether  optical  or  magneto- 
optical.   Unfortunately,  absorption  spectra  involve  both  a 
ground  state  and  an  excited  state,  and  evidence  is  here  pre- 
sented that  the  energy  ordering  of  the  ground  state  MO' s  does 
not  necessarily  determine  the  energy  ordering  of  the  excited 
states,  and  vice-versa.   We  will  return  to  this  after  a  closer 
look  at  the  ground  state. 

The  calculation  shows  the  electronic  structure  built  up 
starting  with  the  "core"  electrons;  Is,  2s,  3s,  2p,  5p  on  Ni 
and  Is,  2s,  on  P.   The  valence  shell  begins  with  the  3d  orbital 
MO's,  so  called  because  Mulliken  population  analysis^  '    shows 
them  to  be  predominantly  metal  3d.   Then  come  all  24  P  2p 
electrons  sandwiched  within  a  narrow  range  of  6,000  cm   .   This 
completes  the  occupied  levels.   The  first  unoccupied  MO  is  the 
Ni  4s,  then  the  3d  2  2,  and  the  4p  last. 

^  ~y 

The  population  analysis  shows  that  the  MO's  are  rather 

completely  segregated  between  the  Ni  and  F  orbitals,  except 

for  the  6a.-,    (3d  2)  MO  which  has  a  substantial  amount  of  P  2p 
Ig  ^  z    ' 

2 
character.   Note  also  that  z   is  the  highest  occupied  3d  orbital 
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and  that  xy  is  lowest.   From  the  orbital  energy  diagram  pre- 
sented here  we  would  expect  the  first  electronic  transition  to 
be  F  2p  to  Nl  4s,  but  In  fact,  as  we  will  show,  this  expectation 
Is  not  fulfilled:   Although  the  orbital  energy  gap  between 
occupied  and  unoccupied  d-orbltal  MO's  Is  over  190,000  cm   , 
the  d-d  bands  are  calculated  to  fall  lowest.  In  the  visible 
and  near  UV  regions  of  the  spectrum. 

The  total  orbital  populations  In  this  computation  are: 
Nl  (core)  18.003^8.13^^0.20^^0.49  ^^^  ^  (core)^-^%^-^\      The 
point  of  Interest  here  Is  that  the  total  occupancy  of  the  Nl  3d 
orbltals  Is  almost  Identical  with  the  free  Ion  configuration, 

o 

3d  .   However,  the  4s  and  4p  orbltals  draw  electrons  from  the 
ligands  so  that  the  computed  net  charge  of  +1.12  on  the  Nl  Is 
substantially  reduced  from  the  formal  Ionic  value  of  +2,  while 
each  F  has  a  computed  net  charge  of  -O.78. 

The  computed  MO's  can  be  Interpreted  to  give  two  types  of 
Information.   The  first  Is  covalency,  which  Is  the  extent  of 
mixing  or  Interaction  of  the  llgand  and  metal  AO's  and  Is 
reflected  In  the  relative  weights  of  metal  and  llgand  coeffi- 
cients In  the  computed  MO's.   The  second  type  of  Information  Is 
the  distortion  of  the  basis  set  AO's  themselves;  I.e.  the  change 
In  the  spatial  extension  of  an  AO  In  the  complex  when  compared 
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to  its  spatial  extension  In  the  free  atom  or  ion.   It  is  un- 
fortunate that  heretofore  all  discrepancies  in  the  properties 
computed  assuming  the  3d  radial  function  of  the  free  Ion  to  be 
a  valid  representation  of  the  3d  AO  in  the  complex  have  been 
attributed  entirely  either  to  expansion  (or  contraction)  of  the 
free  ion  AO,  which  is  the  crystal  field  view,  or  to  covalency 
mixing,  which  is  the  traditional  MO  treatment  of  ESR,  NMR,  and 
optical  data.   In  fact,  both  these  mechanisms  operate  simul- 
taneously and  with  the  use  of  the  double  zeta  representation 
we  can  examine  the  inherent  expansion  or  contraction  of  AO's 
due  to  complex  ion  formation.   This  is  illustrated  in  Table  II. 

The  double  zeta  basis  for  the  F  2p  AO  was  obtained  by 
breaking  up  a  four  term  Gaussian  representation  into  a  long 
range  and  a  short  range  component.   The  expansion  coefficients 
for  atomic  P  and  the  P~  ion  using  the  same  set  of  exponents 
are  given  in  the  first  two  columns  and  beneath  them  is  the  ratio 
in  various  MO's.   Gh  is  the  coefficient  of  the  Gaussian  with 
the  smallest  exponent  so  that  a  smaller  ratio  means  a  more 
diffuse  orbital.   Note  that  the  ratios  for  F  and  F  differ  by 
about  30%. 

The  corresponding  coefficients  in  two  ligand  2p  MO's  are 
shown;  the  tt   non-bonding  Ibp  ,  and  the  slgma  bonding  73-,  „•   The 


ratio  of  Cy'Ch  for  both  these  MO's  is  extremely  close  to  that 

found,  in  the  free  F   ion  calculation. 

2+ 

An  analagous  result  for  the  Ni   3d  AO  is  shown  in  Table 

III.   Here  a  five  term  Gaussian  expansion  was  broken  up  into 
a  set  of  three  and  a  set  of  two.   Here  we  compare  the  ratio 
Cy'Ch  which  does  not  depart  significantly  from  that  computed 
for  the  free  metal  ion  in  its  formal  +2  oxidation  state,  even 
for  the  3d  2  1^0  which  contains  a  significant  amount  of  F~  2p 
character.   Note,  however,  that  the  orbital  deformation,  as 
reflected  in  the  cy'Ch  ratios  is  slightly  different  for  each 
d-orbital  MO.   Thus,  using  a  double  zeta  representation  is 
equivalent  to  optimizing  the  spatial  extension  of  the  basis 

function  separately  for  each  computed  MO. 

2" 
The  excited  electronic  states  of  the  NiFh    ion  are  pre- 
sented in  Table  IV.   The  qualitative  features  of  the  computed 
absorption  energies  are  quite  clear.   The  ordering  of  the  states 
arising  from  d-d  transitions  in  the  virtual  orbital  approxima- 
tion gives  those  originating  from  xy,  and  xz ,  yz  low  in  energy 

2 
and  that  from  z  high  in  energy,  quite  the  opposite  from  that 

expected  from  the  one  electron  orbital  energies  in  the  ground 

state.   The  triplets  are  very  low-lying,  so  low  as  to  fall  in 

the  near-IR  region  of  the  spectrum. 
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2   2 
Since  the  xy,  xz,  yz ,  and  x  -y  MO's  are  overwhelmingly  of 

p- 
metal  character  we  would  expect  that  in  going  from  NiPh    to 

2- 
PdGl|,    only  the  Ni  to  Pd  shift  should  matter,  giving  a  slight 

blue  shift  in  the  computed  energies  of  these  transitions.   On 

2 
the  other  hand,  transitions  originating  from  the  z  MO,  which 

contains  23^  F  2p,  and  considering  the  tremendous  blue  shift 

2~ 
effect  of  the  F   ion,  should  be  red  shifted  in  going  to  PdClj, 

Thus,  the  d-d  bands  in  PdClh    should  be  found  in  the  10-36,000 

-1  2     2   2 

cm   range  with  the  singlet  z  — >x  -y  transition  highest  in 

energy.   These  conclusions  are  not  inconsistent  with  the 

2~ 
observed  absorption  spectrum  of  PdClh   ,  with  the  d-d  bands 

-ifsl 

reported  in  the  range  17-30,000  cm  ^-^'.      However,  a  more  de- 

2~ 
tailed  assignment  of  the  d-d  states  of  PdClh    in  light  of  a 

2- 
calculation  on  NiF^    would  probably  tax  the  powers  of  the 

calculation  beyond  its  present  capabilities. 

Two  high  intensity  bands  have  also  been  observed  in  the 

2  _  _! 

absorption  spectrum  of  the  PdClh    ion  at  36,000  and  ^5>000  cm   , 

with  oscillator  strengths  estimated  at  about  .2  and  .5,  respec- 

(7) 
tively.   Anex^  '    has  measured  the  polarizations  of  these  bands 

and  finds  them  both  to  be  xy  or  in-plane  polarized,  with  a  weak 

z_  polarized  band  hiding  under  the  higher  energy  xy  polarized 

band . 
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The  lowest  energy  symmetry  allowed  transitions  calculated 

2- 

for  NlP^i^   are  shown  in  Table  V  together  with  their  oscillator 

strengths  as  computed  from  the  dipole  velocity  a;  d  dipole  length. 

These  excited  states  are,  of  course,  computed  to  fall  at  very 

2- 
high  energy  in  NiFh    due  to  the  tremendous  blue  shift  effect 

of  F  .   But  most  surprising  ,  the  assignments  of  the  lowest - 

lying  of  these  transitions  are  all  2p-4s  rather  than  the  usually 

assumed  2p-3d.   The  relative  orderings  and  intensities  of  the 

three  computed  2p-4s  bands  are  consistent  with  the  observations 

of  Anex,  although  the  same  can  be  said  about  the  three  2p-3d 

bands  which  are  next  highest  in  energy.   The  use  of  ar  unopti- 

mized  single  zeta  representation  for  the  Ni  4s  orbital  reiiders 

any  conclusions  about  these  assignments  tentative. 

In  Table  VI  we  give  a  detailed  breakdown  of  the  computation 
of  the  excitation  energies.   The  equation  is  given  at  the  top 
and  the  values  of  the  contributing  terms  for  the  three  single 
d-d  transitions  and  the  lowest -lying  charge  transfer  transition 
are  given  in  their  respective  columns,   e.  and  e.  are  the  or- 
bital  energies,  J.  .  and  K.  .  are  the  Coulomb  and  Exchange  Irtegrals, 
respectively. 

The  importance  of  the  electron  repulsion  correction  of  the 
Coulomb  integral  to  the  simple  differerjces  of  ground  state 
orbital  energies  is  immediately  obvious.   Note  that  for  the  d-d 
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bands  the  Coulomb  integral  is  very  similar  to  that  calculated 
from  the  free  ion  metal  3d  orbitals  alone,  given  in  parentheses 
next  to  the  MO  values.   In  contrast  note  the  small  Coulomb 
integral  value  for  the  2p-'4-s  transition.   Thus  although  the 
orbital  energy  differences  for  the  d-d  transitions  are  consider- 
ably larger  than  those  for  the  charge  transfer  bands,  the  enor- 
mous release  in  energy  obtained  by  splitting  up  a  pair  of  d 
orbital  electrons  more  than  compensates  for  this  large  differ- 
ence; ^L*^'  ^^^  enormous  self -repulsion  of  the  3d  electrons  must 
be  considered  as  a  crucial  factor  in  the  gross  ordering  of  the 

low-lying  excited  state  configurations. 

2     2   2 
The  ratio  of  the  Coulomb  integral  for  the  z  — >x  -y 

transition  in  the  molecule  to  its  free  ion  value  is  much  smaller 

than  for  the  other  d-d  transitions.   This  is  a  direct  outcome 

2 
of  the  admixture  of  F  2p  orbitals  in  the  z  MO.   Thus,  although 

2     2   2 
the  orbital  energy  difference  is  smallest  for  the  z  — >  x  -y 

2 
excitation,  reflecting  the  high  energy  of  the  z  MO,  the  value 

for  the  Coulomb  integral  due  to  covalent  mixing  gives  an  excita- 

2     2   2 
tion  energy  for  z  — >  x  -y  highest  of  all  the  d-d  transitions. 

Thus  the  amount  of  covalent  mixing  and  the  accompanying  variable 
reduction  in  the  electron  repulsion  integrals  apparently  deter- 
mine the  finer  details  of  the  ordering  of  the  excited  states  in 


NIF4 


2- 
The  purpose  of  this  computation  has  been  to  demonstrate  to 
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those  Interested  in  the  application  of  MO  theory  to  problems  in 
Inorganic  chemistry  the  feasibility  of  doing  non-empirical  SCF- 
MO  calculations  on  transition  metal  complexes  with  Gaussian 

type  basis  functions.   The  specific  example  discussed,  square 

2- 

planar  NIFk   ,  although  admittedly  a  hypothetical  complex,  was 

chosen  as  a  prototype  to  illustrate  certain  general  characteris- 
tics of  the  calculation.   Since  the  presentation  of  these  results, 

we  have  redone  the  computations  using  completely  optimized  atomic 

2+ 
wavef unctions  for  both  the  Ni   and  F   ions.   The  results  are  in 

qualitative  agreement  with  those  reported  above.   In  addition, 

SCF  calculations  have  been  completed  on  octahedral  NiF/-   as  part 

of  a  study  of  the  interesting  spin  properties  of  this  species 

which  give  rise  to  transferred  hyperfine  interactions,  spin-orbit 

coupling,  and  electronic  g-factors. 
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Figure    II      GROUND   STATE   OF   NIF,,       (D,,^) 


ih' 


fo   Minority 
Composition 


fo   Majority 
Composition 


Orbital  Energy 
(cm"  ) 


4.2 
8.5 


2p 
2p 
2p 


95.8    4p^ 


92.7 
90.1 


3d  2   2 
X  -y 

4s 


+l4l,o4o 
+139,980 
+114,790 


2p 


24 
electrons 


-  17,790 

-  23,360 


22.9 

2p 

77.1 

5.7 

2p 

94.3 

4.4 

2p 

95.6 

3d   2 
2 


3d 
3d 


xz,yz 


Core 


xy 


34 
electrons 


46,530 
50,370 
67,390 
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Table  I  BASIS  SET  FOR  NIF,, 


GTP's     GTO's 


Nl^-^ 

is. 

is^. 

2s, 

2s^ 

3s, 

4s 

6 

12 

2p, 

2V    , 

3P, 

4p 

12 

24 

5d, 

3d^ 

12 

30 

P" 

Is, 
2p, 

2s 
2P^ 

2 

6 

62 

7 

12 
144 
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Table  II   EXPANSION  AND  CONTRACTION  OP  FLUORINE  AO's 


atomic  (P  )     atomic  (P°)     lbp^^(7r"'^* )     Ta-j^  (a^) 


'2u' 


g 


.02652 

.02905 

.17158 

.18408 

.43180 

.49025 

,60428 


.53497 


,44212 


59240 


,35716 
,51830 


cy'c 


.715 


,916 


746 


,689 
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Table  III   EXPANSION  AND  CONTRACTION  OF  NICKEL  AO's 


atomic  (N1   ) 


xy 


xz,yz 


. 07502 
.23864 
.45326 


,44236    .38260    .43072 


.40293 
.17167 


.39003    .35841    .39801 


'3/^4 


1.125 


1.134  1.067  1.082 
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Table  IV   ELECTRONIC  STATES  OF  NIP^^^ 


One  Electron  Energy^  ^^ 

Transition    State    (lOOO  cm  -^)   PdCl,,   (e^^^.J   PtClj,  i^^n^yj 


A 


Ig 


0 


2   2   3 
X2,yz  ->  X  -y   ^E 


g 


2   2     3„ 
xy  -^  X  -y     ^A 


2g 


^  2   2   1^ 
xz,yz  ^  X  -y    E 


g 


2   2     1 
xy  ->  X  -y^    ^A 


2g 


2.22     3-D 
z  ^  X  -y     ^B 


Ig 


2^22     1^ 
z  ^  X  -y      B 


Ig 


^.0 

? 

18.0 

(2) 

11.0 

? 

20.7 

(2) 

19.0 

22.8    (104) 

29.5 

(46) 

19.5 

20.0    (67) 

26.3 

(28) 

26.0 

17.5    (13) 

24.0 

(5) 

^0.0 

29.5    (67) 

36.5 
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Table  V   CHARGE  TRANSFER  STATES  OF  NIF,,^ 


One  Electron 
Transition 

State 

Energy  (cm   ) 

f(V) 

f(r) 

2p  -^  4s 

\^x,y) 

69,400 

.055 

.045 

2p  -^  4s 

\(x,y) 

86,200 

.491 

.263 

2p  ^  4s 

\u(-) 

91,400 

.195 

.106 

2p  -}   5d 

\(x,y) 

100,100 

.081 

.212 

2p  ^  3d 

S,(z) 

101,600 

.000 

.000 

2p  ^  3d 

\(x,y) 

107,000 

.165 

.285 

3d  ^  4p 

V(^) 

108,000 

.209 

.270 
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Table  VI   EXCITATION  ENERGIES  OF  NIF 


•  T.  2- 


4 


^E  (1  ^j)  =  £   -  £.  -  J   +2K 


1  a.-u.  =  219,450  cm 


-1 


Transition   e.-e.(a.y.) 
J   1 


2   2 
xz,yz  ^  X  -y 


V  2   2 
xy  ->  X  -y 


2     2   2 
z  ^  X  -y 


J.  .(a.u.)    K.  .(a.u.)   Energy(cm   ) 


.86741 


.9^495 


.84991 


84913(. 95806)^^^ 

.03455 

19,200 

89496(. 99^29)^^^ 

.01940 

19,500 

73132 (. 9^598 )(^) 

.03264 

40,300 

2p  ^  4s 


,54432      .22783 


.00738     69,400 


(a)   Computed  from  free  Ion  3d  functions, 
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